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Abstract 

High-precision observational data have confirmed with startling evi¬ 
dence that the Universe is currently undergoing a phase of accelerated 
expansion. This phase, one of the most important and challenging cur¬ 
rent problems in cosmology, represents a new imbalance in the governing 
gravitational equations. Historically, physics has addressed such imbal¬ 
ances by either identifying sources that were previously unaccounted for, 
or by altering the gravitational theory. Several candidates, responsible for 
this expansion, have been proposed in the literature, in particular, dark 
energy models and modified gravity models, amongst others. Outstanding 
questions are related to the nature of this so-called “dark energy” that 
is driving this acceleration, and whether it is due to the vacuum energy 
or a dynamical field. On the other hand, the late-time cosmic accelera¬ 
tion may be due to modifications of General Relativity. In this work we 
explore a generalised modified gravity theory, namely f{R, cj), X) gravity, 
where R is the Ricci scalar, 0 is a scalar field, and X is a kinetic term. 
This theory contains a wide range of dark energy and modified gravity 
models. We considered specific models and applications to the late-time 
cosmic acceleration. 
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1 Introduction 


A central theme in Cosmology is the fact that the Universe is currently undergo¬ 
ing an accelerating expansion mm- In this context, during the last two decades 
Cosmology has evolved from being mainly a theoretical area of Physics to be¬ 
come a field supported by observational data. Recent experiments call upon 
state of the art technology in Astronomy and Astrophysics to provide detailed 
information on the contents and history of the Universe, which has led to the 
measuring of the parameters that describe our Universe with increasing preci¬ 
sion. The standard model of cosmology is remarkably successful in accounting 
for the observed features of the Universe. However, there remain a number of 
open questions concerning the foundations of that standard model. In partic¬ 
ular, we lack a fundamental understanding of the mechanisms underlying the 
acceleration of the late-time universe. What is the so-called “dark energy” that 
is driving the cosmic acceleration? Is it a vacuum energy or a dynamical field? 
Is the acceleration due to modifications of Einstein’s theory of General Rela¬ 
tivity (GR)? How is structure formation affected in these alternative scenarios? 
What happens to the universe in the asymptotic future? 

The resolution of these fundamental questions, looking beyond the standard 
theories of gravity and particle physics, is extremely important for theoretical 
cosmology and theoretical physics as a whole. The standard model of cosmology 
has favoured dark energy models as fundamental candidates responsible for 
the cosmic expansion. However, it is clear that these questions involve not 
only gravity, but also particle physics. String theory provides a synthesis of 
these two branches of physics and is widely believed to be moving towards 
a viable quantum gravity theory, given time. One of the key predictions of 
string theory is the existence of extra spatial dimensions. In the brane-world 
scenario, motivated by recent developments in string theory, the observed 4- 
dimensional universe is embedded in a higher-dimensional spacetime [3]. The 
new degrees of freedom belong to the gravitational sector, and can be responsible 
for the late-time cosmic acceleration mm- On the other hand, generalisations of 
the Einstein-Hilbert Lagrangian, including quadratic Lagrangians which involve 
second order curvature invariants have also been extensively explored iniDH 
i sni [n]. An alternative approach was formulated in [121 IE]- Gurvature- 
matter couplings have also been extensively analysed miiiiiiiiiiiiiiHiiiiiiig. 
While these modified theories of gravity offer an alternative explanation to the 
standard cosmological model for the expansion history of the universe EH, it 
offers a paradigm for nature fundamentally distinct from dark energy models 
of cosmic acceleration |22j , even those that perfectly mimic the same expansion 
history. It is also fundamental to understand how one may differentiate these 
modified theories of gravity from dark energy models. 

In this context, a promising way to explain the late-time cosmic acceleration 
is to assume that at large scales Einstein’s theory of GR breaks down, and a more 
general action describes the gravitational field. Thus, one may generalise the 
Einstein-Hilbert action by including “quadratic Lagrangians”, involving second 
order curvature invariants. Some of the physical motivations for these modifica¬ 
tions of gravity were related to the possibility of a more realistic representation 
of the gravitational fields near curvature singularities and to create some first or¬ 
der approximation for the quantum theory of gravitational fields. On the other 
hand, field theories in theoretical physics tends to have Lagrangians quadratic in 
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field strengths, gravity being somewhat different. One may tackle the problem 
using the metric formalism, which consists in varying the action with respect 
to the metric, although other alternative approaches have been considered in 
the literature, namely, the Palatini formalism m, where the metric and the 
connections are treated as separate variables; and the metric-affine formalism, 
where the matter part of the action now depends and is varied with respect to 
the connection [53]. Recently, a novel approach to modified theories of grav¬ 
ity that consists of adding to the Einstein-Hilbert Lagrangian an f{R) term 
constructed a la Palatini ESj. It was shown that the theory can pass the So¬ 
lar System observational constraints even if the scalar field is very light. This 
implies the existence of a long-range scalar field, which is able to modify the 
cosmological and galactic dynamics, but leaves the Solar System unaffected. 
These explicit models are consistent with local tests and lead to the late-time 
cosmic acceleration, and also verify the absence of instabilities in perturbations. 

In this work, we intend to consider a generalised modified gravity theory, 
f{R, (j), X) gravity, where R is the Ricci scalar, (j) a scalar field, and X a kinetic 
term. This theory contains a wide range of known dark energy and modified 
gravity models, for instance f{R) gravity models or Galileons. Here we apply it 
to the late-time cosmic acceleration. In particular, several cosmological solutions 
are studied within the framework of these theories, specifically solutions that can 
provide cosmic acceleration at late times, and even the exact ACDM evolution. 
Reconstructions techniques are implemented in order to obtain the f(R,(j),X) 
of the action given a particular Hubble parameter. This provides a way to 
efficiently check the viability of any gravitational action by just considering a 
particular cosmological evolution and then analysing the gravitational action. 
An action of the form f{R,(j),X) is in fact quite natural as it removes any 
assumptions on the underlying theory of gravity with the exception of being 
second order. We can think of the field (f) as the effective field controlling the 
strength of the gravitational force. One should emphasise that normal matter 
is still coupled minimally to this theory. 

This paper is organised in the following manner: In Section|5J we present the 
general formalism of f{R,(j),X) gravity. In Section [S] we analyse cosmological 
applications, in particular, to the late-time cosmic acceleration, by consider¬ 
ing specific models and reconstructing the corresponding gravitational action. 
Then, we discuss in the conclusions at section Sj about the results obtained and 
possible future considerations of these theories. 

2 Generalised gravity models: formalism 

Consider the following general action 



( 1 ) 


where / = f{R, 4>, X) is a function of the Ricci scalar i?, the scalar field (/>, and a 
kinetic term X = —(I/2)(V0)^. The matter Lagrangian is Lm = Lm{g,fp), we 
note that matter is minimally coupled to the gravitational sector of the theory 
and V' collectively denotes any matter fields. 

The variation of action m with respect to the metric g^^i, gives to the fol- 
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lowing field equation 

= \{f- RF) - 5^. V„ V“F + i/,x , (2) 

while variations with respect to the scalar field </> yield 


V^(/,x0’^) + /.0 = O. (3) 

Here F = df /dR and is the matter energy-momentum tensor defined by 

^(m) _ ^ Lm) / .N 

~ S{9^^) ■ ^ ^ 

Now, taking into account a flat Friedman-Lemaitre-Robertson-Walker (FLRW) 


metric in spherical coordinates 

= —dt^ + a^(t) [dr^ + {dO^ + sin^ 9dip)\ , (5) 

where a{t) is the scale factor, then Eqs. (H)) and (|31) take the following form 

3Fij2 = /^x + i(FR-/)-3idF + «V„,, (6) 

-2FH = fxX + F-HF + K^{p^+p^), (7) 

0 = ^ - f,4 ,, (8) 


where id = d/a is the Hubble parameter and the overdot denotes a derivative 
with respect to cosmological time t. For simplicity we will not take into account 
additional forms of matter. 

In order to derive an equation for the effective equation of state, we rewrite 


the field equations in the following form 

3H^ = (pm + Pde) , (9) 

— 2H = (pm +Pm + PDE + PDe) , (10) 

where the quantities pde and pde are defined by 

k^Pde = f,xX + ^{FR-f)-3HF + 3H\l-F), (11) 

AC%E = F + 2HF -^(FR-f)- (2H + (1 - F) . (12) 

One may easily show that the these effective dark energy components satisfy 
the usual conservation equation, i.e., 

Pde + 37d (pde+Pde) = 0. (13) 


This is expected since the left-hand sides of (©-(Uni) are equivalent to the general 
relativistic field equations containing the Einstein tensor which satisfies the twice 
contracted Bianchi identities. 

Let us define the dark energy equation of state given by wde = Pde/pde, 
which by using Eqs. (IHll- dTUl) takes the form 


wde = — 


2H + 

3id2 - K^p^ 


(14) 


4 




In order to specify the matter content, let us choose the equation of state Pm 
wpm- In turn, we can write the dark energy equation of state as 


wde = — 



(15) 


where 



(16) 


In the following we will consider the above equations and seek solution cor¬ 
responding to an accelerated expansion. 

3 Cosmological applications: Late-time cosmic 
acceleration 

Let us consider some specific forms of the action o and reconstruct particular 
cosmological solutions capable of reproducing the late-time cosmic acceleration. 
We are considering pure de Sitter solutions, also power law solutions and finally 
we are reconstructing an exact ACDM model. Note that this kind of cosmo¬ 
logical solutions have already been considered in the literature within modified 
gravity theories, specially in f{R) gravity and Gauss-Bonnet gravity. In partic¬ 
ular, power-law solutions were found in [26l [27] while exact ACDM model were 
reconstructed in [28| |29l [30]. In order to simplify the calculations, we restrict 
to a matter Lagrangian with a pressure-less fluid, as natural when analysing 
late-time cosmology. 

3.1 Brans-Dicke type models 

Here we consider a Brans-Dicke type action with f{R, X, (p) = j(X, 4i)R, where 
the coupling to the Ricci curvature also includes the kinetic term of the scalar 
field. 



(17) 


Then, by considering an homogeneous and isotropic metric ®, the FLRW equa¬ 
tions (EUTj) are expressed as follows: 



3^{X,cP)H^ 


2^{X,cP)H 


And combining both equations, a differential equation solely of ^{X{t), Pit)) = 
7 (t) is obtained 



(19) 


Then, equation (HID can be solved by assuming a particular Hubble parameter 
and action (HZD can be reconstructed by assuming a particular form of ^{X, (p), 
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which has to satisfy 


l,xX 


+ SiJy — 
R 


( 20 ) 


As an example for illustrating the reconstruction method, here we consider 
the usual form of a scalar field Lagrangian, 


7(X,<(.) = X-F(</)), (21) 

where V{(j)) is a function of the scalar field to be determined. Then, by redefining 
the scalar field to coincide with the cosmic time (j) = t and using (uni and (EU, 
the kinetic term and the scalar potential are fully reconstructed: 

37i?2 _ ^2p^^-3fd4,H 

R ^ 

= X{ci>) - ^{cj^). ( 22 ) 


X(0) 

V{^) 


where we have used the continuity equation pm + 3Hp^ = 0. 

Let us now consider several cosmological solutions. Firstly we assume a pure 
de Sitter solution, 

H = Ho, (23) 

where Hq is a constant. Then, the equation (ITO)) leads to 

+ + = (24) 

which is the equation of a damped harmonic oscillator, whose general solution 
is given by 

-f{t) = e-^°*CcosV2Ho(t-to), (25) 

where {C,to} are integration constants. Then, assuming (I^Tl) . the gravitational 
action is reconstructed by the expressions (1221) . where the kinetic term and the 
scalar potential turn out 


^(0) 

V{<^) 


c 

2 V 2 

c _ 




Ho4> 


4 cos V2Ho((j) — 0o) + sin \/2Ho{(j) — (j)o) 


K. po ^_3 


12HI 


The interesting point of this action is that leads to a constant Hubble parameter 
in the presence of dust matter and a scalar field but without a cosmological 
constant, contrary to standard General Relativity. 

Another important class of solutions in cosmology is the power-law type 
expansion, whose Hubble parameter is given as follows 

71 

H=- ^ a{t) = aor. (27) 


By applying the same procedure as above, the equation (HU) turns out to be 


dt^ 


+ 2nt—^^ + n{3n — 2)^{t) = 0 , 


(28) 
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which is an Cauchy-Euler equation whose general solution is given by 


7(t) = CitP+‘> + C2tP-^ , 


( 29 ) 


where Ci and C 2 are two constants of integration, and we defined 


P = 


1 - 2 n 


q = 


\/l + 4n — 8 n^ 


(30) 


By assuming (EU as above, the kinetic term and the scalar potential are 
obtained 


VW 


1 

4(2n- 1) 

1 

4(2n- 1) 


Ci(l + + C2{1 - 2q)r-^ - , 

Cl(5 - 8 n + 2q)^+‘i + C 2(5 - 8 n - 2g)0P-9 - 


Hence, the gravitational action capable of reproducing the power-law solutions 
(Ell) is obtained. 

Finally, let us consider the case of exact ACDM and explore whether a 
gravitational action of the type (EZD is capable of reproducing such kind of 
expansion. In order to simplify the calculations, here we use the redshift 1 + z = 
i as the independent variable instead of the cosmic time t. The cosmological 
evolution of ACDM model can be expressed by the following Hubble parameter 

H{z) = + z)3 + 1 - , (32) 


where ilm = 3^^°/^2 and Hq is the Hubble parameter evaluated at z = 0. Then, 
the equation m can be rewritten in terms of the redshift as, 

(1 + z)"hV' + [(1 + z fHH' - (1 + z)H^] 7' + [m^ - 2(1 + z)HH'] 7 = 0 , 

(33) 

where the primes denote derivatives with respect to the redshift. Hence, the 
equation (I!I1I)) can be solved as 'y{X,(j)) = 7(2). Particularly, for the ACDM 
model, the equation dMl) leads to 

2 ( 1 + 2 :)^ + 3z “h 3) + 1 ] 7'^^(z)H-(1+z) “t- 3z^ H- 3z + 3) — 2j ^^(z)+ 6 (l—i277j)'y(z) = 0 . 

(34) 

Unfortunately, this equation can not be solved exactly but numerical re¬ 
sources are required. In order to illustrate a particular case, we solve the equa¬ 
tion (IMl) by assuming Dm = 0.31, provided by the last Planck data EH- The 
kinetic term and the scalar potential are depicted in Fig. [TJ Note that we have 
assumed the initial conditions "f{z = 10) = 1 and 7 '(z = 10) = 0, so that GR 
is recovered at large redshifts. As expected, the dynamics of the scalar field 
become important at small redshifts while remains virtually constant at large 
ones. 


3.2 Minimally coupled field model 

Let us now consider another different type of the gravitational action CD, where 
the scalar field is minimally coupled. We choose f{R,X,cj)) = a{R) +'-f{X,(j)) 
so that 

S = J dx‘^^/^[a{R)+j{X,(j)) + L,^] . (35) 
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Figure 1: The kinetic term (left panel) and scalar potential (right panel) for the 
ACDM model (15^ with Vim = 0.31. 


This action can be seen as an /(i?)-like action with a k-essence-like fluid in the 
presence of f{R) gravity. Then, the FLRW equations ([HllTl) turn out 

+^{a^nR-a-'y) - 3Ha^R +K^pm, 

-2a^RH = ^^xX + a^R - H6 i^r + K^pm, (36) 

which can be combined, so that an equation for ^(X, (j)) = jlyt) is obtained, 
j{t) = -a + (i? - - 4i/) a,R - (diJi? + 2R^ a,RR - 2R^aRRR . (37) 


Then, by specifying a{R) and providing a particular Hubble evolution H{t), 
the function ^{t) is determined. The complete function j^X, 4>) is provided by 
assuming a particular form of the scalar field Lagrangian that satisfy 


1 


1 


1 


l,xX = -a + 3i7 - -R]a^R + SHa^R - k Pm + xl ■ 


(38) 


In order to reconstruct some particular solutions, we are considering the follow¬ 
ing a{R), 

a{R)=R + aR^, (39) 

where a is a constant free parameter with the appropriate dimensions. This 
gravitational action is a very well known one in the literature as it is capable of 
reproducing inflation [32] ■ As above, we start by considering a pure de Sitter 
solution (1331) . H = Hq. In such case, the expression of j{t) in (1371) yields 

7 = -6i7o" ■ (40) 

We can assume again 'y{X, </>) = X — V {(j)) for illustrating the full reconstruction 
of the action dSSI), and redefine the scalar field as (j) = t. Then, the kinetic term 
X{(p) and the scalar potential V{(j)) turn out 


X{(1)) = Xoe-^^°'^ 


v{4>) = Xoe-3^°^ + ml , 


(41) 






where Xq = —k^pq. 

In the same way as above, power-law solutions can be easily reconstructed. 
By assuming H{t) = j, the expression (1571) leads to 


7 (t) = 2n 


18o(4 — lln -I- 6n^) -I- (2 — 3n)t^ 


(42) 


while the kinetic term and the scalar potential yield 
A-3n , „a{36-72n) + (j)‘^ 


= Xod'”'"” + 2n 

V{^) = Xor'^^ + 


-3n I —36au(2-I-n(—7-I-6n))-I-2ri(—1-I-3n)(/>^ 


, (43) 


where we have again assumed 7 (X, (j)) = X — V{(f)) for illustrative purposes. 

Finally, let us reconstruct the gravitational action for ACDM model. In such 
a case, the Hubble parameter can be expressed as follows 


H{t) = 





(44) 


where A and to are constants. Then, by using the equation (I37L the expression 
for 7 (t) yields 


lit) 


—2A — 9a A^ sech"^ 


to) 


(45) 


While assuming 7 (X, (f>) = X — V{(j)), the kinetic term and the scalar potential 
result 


X{<j)) 

V{cj)) 


—9aA^ sech'* 


1 

2 


y/^{(j) - to) , 


2A. 


(46) 


Hence, the scalar potential reduces to a cosmological constant, while the dynam¬ 
ics of the scalar field compensates the extra terms in the gravitational sector 

dMl)- 

Another interesting set of cosmological models are the so-called scaling so¬ 
lutions, which are characterized by 

— =C', (47) 

Pm 

where (7 is a constant, and p is the energy density of an additional perfect fluid. 
We can rewrite the first FLRW equation (1361) as follows 
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P + Pm ) 


where 



j,xX + -{a^nR - a - ^) - 3Ha^R , 


which leads to the constraint equation: 


(48) 

(49) 


1 


l,xX + - {u^rR - a - 7 ) - 3Ha^R 


= 0 . 


(50) 
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By considering the Hilbert-Einstein action, a{R) = R, together with a kind of 
k-essence field that satisfies pjPm = C, the equation (HSl) can be easily solved 
by using the continuity equation leading to 


while the equation (1501) leads to 


7.xX--7 = C^ 


(51) 


(52) 


And the second FLRW equation (IM)) provides an additional condition over the 
scalar field sector: 




4 — 3k^po 

’ 


(53) 


which together with dSH) allows to reconstruct the scalar field Lagrangian once 
a particular form of 7 is assumed. In order to illustrate the procedure, let us 
consider the usual quintessence field 7 = — V{(/)). Then, the kinetic term 

and the scalar potential yield 


A(0) 

VW 


4 — 3k^po 

3^2 ’ 

k'^Po{6C + 1) — 4 
302 


(54) 


where once again we have redefined the scalar field as the cosmic time, (j) = t. 
Note that the potential (1501) is very well known as one of the possible scalar 
potentials that provides scaling solutions within quintessence models, others 
are the exponential ones V{(j)) oc e'^, see Refs. [331134] . 


4 Conclusions 

In this work, we explored a generalised modified gravity theory, namely, f{R, 4>,X) 
gravity, where R is the Ricci scalar, 0 a scalar field, and X a kinetic term. This 
theory contains a wide range of dark energy and modified gravity models, par¬ 
ticularly higher derivative models, both in the gravitational sector and also in 
the matter one. Some models within this action include f{R) modified gravity 
and Galileons, which have been widely studied in the literature. We have then 
considered some specific Hubble parameters, specifically those which reproduce 
late-time cosmic acceleration and a realistic cosmological evolution, and the 
corresponding /(i?,0. A) action has been reconstructed. 

Using the reconstruction technique studied here, the gravitational action can 
be easily obtained once some restrictions are assumed on the kind of function 
/(i?, 0, A). We note that this action carries several extra degrees of freedom, 
so a particular cosmological solution does not provide a unique action. A de¬ 
generacy problem as in every dark energy model. However we suggest some 
simple actions that reproduce late-time acceleration and even the exact ACDM 
behaviour. Future works should be devoted to studying the viability of the 
different actions within /(R, 0, A) gravity framework. Besides the cosmologi¬ 
cal evolution, one should study cosmological perturbation theory with the view 
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of understanding the growth of structure. Another interesting point may arise 
when analysing screening-like mechanisms, similar to the chameleon mechanism 
in this framework. Note also that the cases of the gravitational action m 
studied here, stand for a minority of the wide class of models included in the 
action (O, some interesting ones to be considered in future papers may be of 
the type fi{R)f 2 {X,<j)). Taking into account various observational facts about 
cosmology, it is hoped that one can restrict the form of the action f{R,(j),X) 
and provide a clearer way for explaining the late-time acceleration and even its 
unification with the inflationary paradigm. 
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